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Abstract 

We show that Zhang's sandpile model (N, [a,b]) on N sites and 
with uniform additions on [a, b] has a unique stationary measure for 
all < a < b < 1. This generalizes earlier results of [6] where this was 
shown in some special cases. 

We define the infinite volume Zhang's sandpile model in dimension 
d > 1, in which topplings occur according to a Markov toppling pro- 
cess, and we study the stabilizability of initial configurations chosen 
according to some measure p. We show that for a stationary ergodic 
measure p with density p, for all p < i, p is stabilizable; for all p > 1, 
p is not stabilizable; for ^ < p < 1, when p is near to ^ or 1, both 
possibilities can occur. 

1 Introduction 



Zhang's sandpile model [M] is a variant of the more common abelian sandpile 
model [3], which was introduced in [T] as a toy model to study self-organized 
criticality. We define the model more precisely in the next section, but we 
start here with an informal discussion. 

Zhang's model differs from the abelian sandpile model on a finite grid 
A in the following respects: The configuration space is [0, 1) A , rather than 
{0, 1, . . . , 2d — 1} A . The model evolves, like the abelian sandpile model, in 
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discrete time through additions and subsequent stabilization through top- 
plings of unstable sites. However, in Zhang's model, an addition consists 
of a continuous amount, uniformly distributed on [a, b] C [0, 1], rather than 
one 'grain'. Furthermore, in a Zhang toppling of an unstable site, the entire 
height of this site is distributed equally among the neighbors, whereas in the 
abelian sandpile model one grain moves to each neighbor irrespective of the 
height of the toppling site. 

Since the result of a toppling depends on the height of the toppling site, 
Zhang's model is not abelian. This means that 'stabilization through top- 
plings' is not immediately well-defined. However, as pointed out in [6] , when 
we work on the line, and as long as there are no two neighbouring unstable 
sites, topplings are abelian. When the initial configuration is stable, we will 
only encounter realizations with no two neighbouring unstable sites, and we 
have - a fortiori - that the model is abelian. 

In [6], the following main results, in dimension 1, were obtained. Unique- 
ness of the stationary measure was proved for a number of special cases: 
(1) a > ~; (2) N = 1, and (3) [a, b] = [0,1]. For the model on one site 
with a = 0, an explicit expression for the stationary height distribution was 
obtained. Furthermore, the existence of so called 'quasi- units' was proved 
for a > 1/2, that is, in the limit of the number of sites to infinity, the one- 
dimensional marginal of the stationary distribution concentrates on a single 
value ^p. 

In the first part of the present paper, we prove, in dimension 1, uniqueness 
of the stationary measure for the general model, via a coupling which is much 
more complicated than the one used in [6] for the special case a > 1/2. 

In Section HI we study an infinite- volume version of Zhang's model in any 
dimension. A similar infinite-volume version of the abelian sandpile model 
has been studied in (TUJ E] and we will in fact also use some of the ideas in 
that paper. 

For Zhang's infinite- volume model in dimension > 1, we start from a 
random initial configuration in [0, oo) z , and evolve it in time by Zhang 
topplings of unstable sites. We are interested in whether or not there exists 
a limiting stable configuration. Since Zhang topplings are not abelian, for 
a given configuration rj e [0, oo) z , for some sequence of topplings it may 
converge to a stable configuration but for others, it may not. Moreover we 
do not expect the final configuration - if there is any - to be unique. Therefore, 
we choose a random order of topplings as follows. To every site we attach an 
independent rate 1 Poisson clock, and when the clock rings, we topple this 
site if it is unstable; if it is stable we do nothing. For obvious reasons we call 
this the Markov toppling process. 

We show that if we choose the initial configuration according to a sta- 
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tionary ergodic measure p with density p, then for all p < |, p is stabilizable, 
that is, the configuration converges to a final stable configuration. For all 
p > 1, p is not stabilizable. For | < p < 1, when p is near to \ or 1, both 
possibilities can occur. 



2 Model definition and notation 

In this section, we discuss Zhang's sandpile on N sites, labelled 1, 2, . . . , N. 
We denoted by X N = [0, oo)^ the set of all possible configurations in Zhang's 
sandpile model. We will use symbols rj and £ to denote a configuration. We 
denote the value of a configuration rj at site x by i] x , and refer to this value 
as the height, mass or energy at site x. We introduce a labelling of sites 
according to their height, as follows. 

Definition 2.1. For rj G AV, we call site x 

empty if r) x = 0, 

anomalous if rj x G (0, |), 

/uZZ if r\ x £ [|, 1), 

unstable if r\ x > 1. 

A site x is stable for 77 if < 77^ < 1, and hence all the empty, anomalous 
and full sites are stable. A configuration rj is called stable if all sites are 
stable, otherwise r/ is unstable. Qn = [0, 1)^ denotes the set of all stable 
configurations. 

By T x {rj) we denote the (Zhang) toppling operator for site x, acting on rj 
and which is defined as follows. 

Definition 2.2. For all rj G Xn such that r\ x > 1, we define 

{0 ifx = y, 

Vy + \f)x if \y-x\ = 1, 
rj y otherwise. 

For all i] such that r\ x < l,T x (r)) = r], for all x. 

In other words, the toppling operator only changes 77 if site x is unstable; 
in that case, it divides its energy in equal portions among its neighbors. We 
say in that case that site x topples. If a boundary site topples, then half of its 
energy disappears from the system. Every configuration in Xn can stabilize, 
that is, reach a final configuration in Qn, through finitely many topplings of 
unstable sites, since energy is dissipated at the boundary. 
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We define the (N, [a,b]) model as a discrete time Markov process with 
state space f2jy, as follows. The process starts at time from configuration 
77(0) = 77. For every t = 1,2,..., the configuration i](t) is obtained from 
T](t — 1) as follows: a random amount of energy U(t), uniformly distributed 
on [a, b], is added to a uniformly chosen random site X(t) G {1, ...,jV}, 
hence P(X(t) = j) = 1/N for all j = 1, . . . , N. The random variables U(t) 
and X(t) are independent of each other and of the past of the process. We 
stabilize the resulting configuration through topplings (if it is already stable, 
then we do not change it), to obtain the new configuration rj(t). By ~E V and 
F' n , we denote expectation resp. probability with respect to this process. 

3 Uniqueness of the stationary distribution 

In Zhang's sandpile model, it is not obvious that the stationary distribution 
is unique, since the state space is uncountable. For the three cases: (1) 
N = 1, (2) a > \ and N > 2, (3) a = 0, b = 1 and N > 2, it is shown in [BJ 
that the model has a unique stationary distribution, and in addition, in case 
(2) and (3), for every initial distribution u, the measure at time t, denoted 
by Vf ,b ' N , converges in total variation to the stationary distribution. In the 
case N = 1, there are values of a and b where we only have time- average 
total variation convergence, see Theorem 4.1 of [Bj. 

In all these cases (except when N = 1) the proof consisted of constructing 
a coupling of two copies of Zhang's model with arbitrary initial configura- 
tions, in such a way that after some (random) time, the two coupled processes 
are identical. We will call such a coupling 'successful', as in [1 3J . Each cou- 
pling was very specific for the case considered. In the proof for the case a > 
and N > 2, explicit use is being made of the fact that an addition to a full 
site always causes a toppling. The proof given for the case a = and 6=1 
and N > 1 can be generalized to other values of b, but a = is necessary, 
since in the coupling we use that additions can be arbitrarily small. In the 
special case N — 1, the model is a renewal process, and the proof relies on 
that. Here is the main result of this section; note that only the case a = b is 
not included. 

Theorem 3.1. For every < a < b < 1, and N > 2, Zhang's sandpile model 
(N,[a,b\) has a unique stationary distribution which we denote by fi a,b,N . 
Moreover, for every initial distribution on [0, 1)^, the distribution of the 
process at time t converges in total variation to fi a,b ' N , as t — > 00. 

We introduce some notation. Denote 77, £ G as the initial configura- 
tions, and rj(t), £(£) as two independent copies of the processes, starting from 
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r\ and £ respectively. The independent additions at time t for the two pro- 
cesses starting from 77, £ are [/''(i) and U^(t), addition sites are X v {t),X^(t) 
respectively. Often, we will use 'hat '-versions of the various quantities to 
denote a coupling between two processes. So, for instance, f)(t),£(t) denote 
coupled processes (to be made precise below) with initial configurations 77 
and £ respectively. By X v {t) and X^(t) we denote the addition sites at time 
t in the coupling, and by U v (t) and U^(t) the addition amounts at time t. 

In this section, we will encounter configurations that are such that they 
are empty at some site x, 1 < x < N, and full at all the other sites. We denote 
the set of such configurations £ x . By £&, we denote the set of configurations 
that have only one empty boundary site, and are full at all other sites, that 
is, £ b = Si U £ N . 

In order to give the proof of Theorem 13. 1\ we need the following three 
preliminary results, the proof of which will be given in Sections 13.11 13.21 and 
13.31 respectively. 

Lemma 3.2. For all r], £ in Qn, v(t) and £(£) are a.s. simultaneously in Sb 
infinitely often. 

Lemma 3.3. Let 77 and £ be two configurations in Sn and let, for all e > 0, 

Consider couplings (r/(t),£(t)) of the process starting at 77 and £ respectively. 
Let, in such a coupling, T be the first time t with the property that 

max .1 fj x {t) -i x (t) |< e (3.1) 

l<x<N 

and 

j)(t)e£ N ,i(t)e£ N . (3.2) 

There exists a coupling such that the event T < t e has probability at least 
(2N)~ t % uniformly in rj and£. 



Lemma 3.4. Let 



b — a 



e a ,b,N - 1 ctjN 1 



6 + i6IL^ 1 (l + 2 iV - 2 -') 

Consider couplings (r/(t),£(t)) of the process starting atrj and £ respectively, 
with the property that 

max I 773. — £ x \< e aAN . (3.3) 

l<x<N 
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Let T' be the first time t with the property that r)(t) = £(t). Then there exists 
a coupling such that the event T' < (N — 1) f^rg] has probability bounded 
below by a positive constant that depends only on a, b and N. 

We now present the coupling that constitutes the proof of Theorem 13.11 
making use of the results stated in Lemma 13.2} Lemma 13.31 and Lemma 13.41 

Proof of Theorem \3.1[ Take two probability distributions ui, v 2 on Q^, and 
choose r] and £ according to v x , u 2 respectively. We present a successful 
coupling {r)(t), £(£)}, with r/(0) = r\ and £(0) = £. If we assume that both 
v\ and v 2 are stationary, then the existence of the coupling shows that v\ = 
i>2 = v. If we take vi = v and z/ 2 arbitrary, then the existence of the coupling 
shows that any initial distribution z/ 2 converges in total variation to v. 
The coupling consists of three steps, and is described as follows. 

• step 1. We evolve the two processes independently until they encounter 
a configuration in Sb simultaneously. From Lemma 13.21 we know this 
happens a.s. By symmetry, we can assume that both configurations 
are in Sn as soon as both processes have reached a configuration in 
From that moment on, we proceed to 

• step 2. We use the coupling as described in the proof of Lemma 
1331 That amounts to choosing X*(t) = X*(t) = Xi(t), and U^(t) = 
tr>(t) = U^{t). As the proof of Lemma shows, if U^{t) and X*(t) 
satisfy certain requirements for at most t e time steps, then we have that 
(13.11) and (13.21) occur, with e = e^^. If during this step, at any time 
step either U v (t) or X v (t) does not satisfy the requirements, then we 
return to step 1. But once we have ( 13. ip and (13.21) (which, by Lemma 
13.31 has positive probability), then we proceed to 

• step 3. Here, we use the coupling as described in the proof of Lemma 
EH Again, we choose = X r '(t) = X"(t) and t/"(t) = lF(t), but 
the dependence of U^(t) on U v (t) is more complicated; the details can 
be found in the proof of Lemma |3~4"1 As the proof of Lemma l3~4l shows, 
if U v (t) and X^it) satisfy certain requirements for at most (iV— 1) [^4^] 

time steps, then we have that fj(t) = occurs, and from that moment 
on the two processes evolve identically. By Lemma 13. 4[ this event has 
positive probability. If during this step, at any time step, either U n {t) 
or X v (t) does not satisfy the requirements, we return to step 1. 

In the coupling, we keep returning to step 1 until step 2 and subsequently 
step 3 are successfully completed, after which we have that fj(t) = Since 
each step is successfully completed with uniform positive probability, we a.s. 
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need only finitely many attempts. Therefore we achieve rj(t) = in finite 
time, so that the coupling is successful. □ 

Now, we will proceed to give the proof of Lemma 13.21 Lemma 13.31 and 
Lemma 13.41 

3.1 Proof of Lemma [3721 

In this section, we show that starting two independent processes from any 
two configurations r) and £, the two processes will a.s. be in Sb simultaneously 
infinitely often. The proof will be realized in two steps. 

Lemma 3.5. Let 77 be a configuration in Qn- The process starting from t] vis- 
its £n within (N + 1) \-^\ time steps, with probability at least (^)( N+1 ^a+b\ 

Proof. We prove this by giving an explicit event realizing this, that has the 
mentioned probability. In this step, we always make heavy additions to site 
N, that is, additions with value at least (a + b)/2. 

First, starting from configuration rj, we make heavy additions to site 1 
until site 1 becomes unstable. Then an avalanche occurs and a new con- 
figuration with at least one empty site is reached. The leftmost empty site 
denoted by r%. If r\ = N we are done. If t\ 7^ N, then it is easy to check 
that site r% + 1 is full. The total number of additions needed for this step is 
at most 2r^-rl. 

I a+b I 

Then, if r\ 7^ N, we continue by making heavy additions to site r\ + 1 
until site r\ + 1 becomes unstable. Then an avalanche starts from site ri + 1. 
During this avalanche, sites 1 to r\ — 1 are not affected, site T\ becomes full 
and we again reach a new configuration with at least one empty site, the 
leftmost of which is denoted by r 2 . If r 2 = N we are done. If not, note that 
r 2 > r i + 1 and that all sites 1, r 2 — 1 and r 2 + 1 are full. At most 
heavy additions are needed for this step. 

If r 2 ^ N, we repeat this last procedure. After each avalanche, the 
leftmost empty site moves at least one site to the right, and hence, after the 
first step we need at most N — 1 further steps. 

Hence, the total number of heavy additions needed for the above steps 
is bounded above by (N + Every time step, with probability a 

fixed site is chosen and with probability |, an addition is a heavy addition. 
Therefore, the probability of this event is at least (2N)^ ( - N+1 ^~\ □ 

Lemma 3.6. Let £(0) G £&, then £(1) e Sb with probability at least j?. 
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Proof. Again, we give an explicit possibility with probability -k. Starting in 
£ G £b, we make one addition to the site next to the empty boundary site. 
If this site does not topple, then of course £(1) is still in Sb. But if it does 
topple, then every full site will topple once, after which all sites will be full 
except for the opposite (previously full) boundary site. In other words, then 
£(1) is also in S^. The probability that the addition site is the site next to the 
empty boundary site, is j?. Then £(1) G with probability at least 4. □ 

Proof of Lemma \3.S\ From Lemma 13.51 it follows that the process starting 
from £ is in Sb infinitely often. Let be the first time that the process is in 
Sb, and define 

T b v = mm{t:t>tlr](t) G S b }. 

By the same lemma, the probability that < T b v — < (N + 1) [^tj] is at 

least (2N)-( N+1 ^^\ 

Repeatedly applying Lemma l3~o1 gives that the event that + 1) G 
£b,£(tl + 2) G £b, ...,£(T^) G £b, occurs with probability bounded below by 

We have showed that when G £b, within at most (N + 1) [^tj] time 
steps, the two processes are in £b simultaneously with probability at least 
(l/W Y N+1 ^~^ . Combining this with the fact that the process starting from 
£ is in £ b infinitely often, we conclude the two processes are in £ b simultane- 
ously infinitely often. □ 



3.2 Proof of Lemma [3731 

In this part, we couple two processes starting from 77, £ G £jy. The coupling 
consists of choosing the addition amounts and sites equal at each time step. 
For this coupling, we present an event that has probability (2iV)~* £ , with 
t e = 2\-^] ■ [log jj-r^i) (ff)l' an d which is such that if it occurs, then 
(13.11) and ( 13.2ft are satisfied. 

The event we need is that for t t time steps, 

1. All additions are heavy. 

2. The additions occur to site iV until site iV becomes unstable, then to 
site 1 until site 1 becomes unstable, then to site iV again, etcetera. 

The probability for an addition to be heavy is ~ and the probability for the 
addition to occur to a fixed site is -h. Therefore, the probability of this event 
is (2iV)-* e . 
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Now we show that if this event occurs, then (13. ip and (13. 2p are satisfied. 
Let U(t) be the addition amount at time t. Define a series of stopping times 
{rk}k>o by 

r = 0, r k := min{t > r fc _i : ^ U(t) > 1}, for k > 1, (3.4) 

t=T fc _l+l 

and denote 

Tk 

s k = Yl ^(*)- ( 3 - 5 ) 

The times (/c > 0) are such that in both configurations, only at these 
times an avalanche occurs. Indeed, for the first avalanche this is clear because 
we only added to site N, which was empty before we started adding. But 
whenever an avalanche starts at a boundary site, and all other sites are full, 
then every site topples exactly once and after the avalanche, the opposite 
boundary site is empty. Thus the argument applies to all avalanches. 

Since we make only heavy additions, 

2 

U - Tfc-i < [ — — 1, for all k. (3.6) 
a + b 

After the k-th avalanche, the height %{rk) is a linear combination of 
Si, Sk and 7/i, r/Tv-i, which we write as 

k N 

Vy{n) = ^2 A i y {k)Si + Y B my{k)rim, for 1 < y < N, (3.7) 

1=1 m=l 

and a similar expression for ^(t^). From Proposition 3.7 of [6], we have that 
B my {k) < {1 - 2~^) max B mx {k-1). 

X 

By induction, we find 

B my {k) < (1 - 2-^) k 

and hence 

max | f) y (T k ) - iy(T k ) I < N(l- 2-r¥l) fc max | r] x - & \ 
iSyS" i<x<n 

where we use the fact that i],t; e En implies m&xi< x <N \ t\ x — £ x \< \. 
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For each e > 0, choose k e = 2 flog (1 _ 2 _ r ^i Then f (1-2~^1)^ < 
e, so that 

max | 7) x (TfcJ - LinJ |< e, 

l<a:<iV 

and moreover, an even number of avalanches occurred, which means that at 
time Tfc e , both processes are in Sn- By (13.61) . Tk c <t e = k e \-^\. Thus, is 
a random time T as in the statement of Lemma 13.31 □ 

3.3 Proof of Lemma 13.41 

As in the proof of Lemma 13.31 we will describe the coupling, along with an 
event that has probability bounded below by a constant that only depends 
on a, b and N, and is such that if it occurs, then within (JV — 1) TsTfel ^ me 
steps, fj(t) = £(£). First we explain the idea behind the coupling and this 
event, after that we will work out the mathematical details. 

The idea is that in both processes we add the same amount, only to site 1, 
until an avalanche is about to occur. We then add slightly different amounts 
while still ensuring that an avalanche occurs in both processes. After the first 
avalanche, all sites contain a linear combination of the energy before the last 
addition, plus a nonzero amount of the last addition. We choose the difference 
D\ between the additions that cause the first avalanche, such that site N will 
have the same energy in both processes after the first avalanche, where \D\\ 
is bounded above by a value that only depends on a, b and N. Sites N — 1 
will become empty, and the differences between the two new configurations 
on all other sites will be larger than those before this avalanche, but can be 
controlled. 

When we keep adding to site 1, in the next avalanche only the sites 
1, . . . , N — 2 will topple. We choose the addition amounts such that after the 
second avalanche, sites N — 1 will have the same energy. Since site N did 
not change in this avalanche, we now have equality on two sites. After the 
second avalanche, site N — 2 is empty, and the configurations are still more 
different on all other sites, but the difference can again be controlled. 

We keep adding to site 1 until, after a total of iV — 1 avalanches, the 
configurations are equal on all sites. After each avalanche, we have equality 
on one more site, and the difference increases on the nonequal sites. We 
deal with this increasing difference by controlling the maximal difference 
between the corresponding sites of the two starting configurations by the 
constant e. a ,b,Ni so that we can choose each addition of both sequences from 
a nonempty interval in [a, b\. The whole event takes place in finite time, and 
will therefore have positive probability. 
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Proof of Lemma \3.4\ The coupling is as follows. We choose X v (t) = X n (t), 
and £>(t) = U^{t). We choose the addition sites Xi(t) and X^(t) equal, 
and the addition amounts U v (t) and U*(t) either equal, or not equal but 
dependent. In the last case, U^(t) is always of the form a + (U v (t) + D — a) 
mod {b — a), where D does not depend on U v (t). The reader can check that, 
for any D, U^(t) is then uniformly distributed on [a, b}. 

The event we need is as follows. First, all additions are heavy. For the 
duration of N — 1 avalanches, which is at most (N — 1) ["^rg] time steps, the 
additions to i] occur to site 1, and the amount is for every time step in a 
certain subinterval of [a, b], to be specified next. 

We denote ^ = a" and recursively define 

e k+1 = (l + 2 N ~ k ~ 2 )e k , 

with ex = e a ,b.N- Between the (k — l)-st and k-th avalanche, the interval is 
[a", a" + 2efc], and at the time where the k-th avalanche occurs, the interval 
is a subinterval of [a 1 , b] = [a" + 3e k , b] of length see below. 

The probability of at most (N — 1) additions occurring to site 1, 

is bounded from below by jV~^"+&""- Ar ~ 1 ' ) . Since e k is increasing in k, the 
probability of all addition amounts occurring in the appropriate intervals, is 
bounded below by (2ei) (JV - 1)( W -1) ■ - The probability of 

the event is bounded below by the product of these two bounds. 

Now we define the coupling such that if this event occurs, then after N — l 
avalanches, we have that f)(t) = £(£). In the remainder, we suppose that the 
above event occurs. 

We start with discussing the time steps until the first avalanche. Suppose, 
without loss of generality, that 771 > £1. We make equal additions in [a", a" + 
2ei], until the first moment t such that f]\{t) > 1 — a" — 2t\. We then 
know that > 1 — a" — 3ei. If we now choose the last addition for both 
configurations in [a', b] = [a" + 3ei, b], then both will topple. 

Define 

N-l 

D 1 :=Y,^\Vy-Q- (3-8) 

y=X 

Let 7~i be the time at which the first avalanche occurs. Then we choose for 
all t < n, fpi(t) = U*(t), and L^(n) - U^n) = D v Since |£>i| < 2 N ~ l e 1 < 
when U*(ti) E ^] (the middle half of [a', b]) 

a' < U v (n) + D x <b. (3.9) 

So, if £>?(ti) is uniformly distributed on [3s£t6, 2^36] j then U^n) = U^n) + 
D x is uniformly distributed on [2^6 + D u + Dx\ C [a', 6]. 
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Let Ri = YH=\ U v (t). Then at time Ti, for 1 < x < N — 2 we have 

Un) = ^t(vi + Ri + u ri (n)) + -V2 + --- + \n x+1 , 

and 

Vn-i(ti) = 0;57at(ti) = fj N _ 2 (Ti), 
and a similar expression for £ x (tl)- It follows that 

iV-1 

^(ri) - Un) = -2^D 1 + £ ^(Vy -Q = 

y=l 

which means that the two coupled processes at time T\ are equal at site N. 

After this first avalanche, the differences on sites 1, . . . , N — 3 have been 
increased. Ignoring the fact that sites N — 2 happen to be equal (to simplify 
the discussion), we calculate 



max I 77 x (ti) — £ x (ti) |< max <( — — I D\ I + max 

\<x<N 1 ' ^ ' S V ' l<x<N I 2 X+1 ' ' 1<^<JV 



x+1 x \ 



f2 iV - 1 1 I 

< max < - ■ - + (1 — — ) )■ max \ri v — £ v \ 

~ l<x<N \ 2 X+1 V 2 X ' J l<x<N 1 W SJ/ 1 

<(l + 2^ 3 ) max 

l<a:<iV 

<(l + 2^ 3 ) ei = e 2 . (3.10) 

We wish to iterate the above procedure for the next N — 2 avalanches. We 
number the avalanches 1, . . . , N — 1, and define as the time at which the k- 
th avalanche occurs. As explained for the case k — 1, we choose all additions 
equal, except at times r^, where we choose U^ijk) — U^ijk) = Dk, with 

N-k 

D k = Y,^ l [%^i)-Ur k -i)] 

y=l 

and 

\D k \ < 2 N ~ k niax ^(rfc-i) -| x (r fc _i)|. 

l<x<N 

The maximal difference between corresponding sites in the two resulting 
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configurations has the following bound: 



max - \rfx(r k ) ~ 

Kx<N 



x+1 ^ 




Z l<x<N 



max |?7 x (T fe _i) - ^ x (r fc _i' 

1<I<JV 



(3.11) 

Hence for all k, \D k \ is bounded from above by e k+1 2 N ~ k , where e k+ i only 
depends on e k and N. With induction, we find, 



|D Jfc |<2 JV -*nf=i(l + 2 JV - l - 2 )ci:=d*. 



We choose e\ = e a ,b,N such that D N _i < As the upper bound d k is 

increasing in k, we get D k < for all k — 1, AT — 1. □ 

Remark It follows from our proof that the convergence to the stationary dis- 
tribution goes in fact exponentially fast. Indeed, every step of the coupling 
is such that a certain good event occurs with a certain minimal probability 
within a bounded number of steps. Hence, there exists a probabiliy q > 
and a number M > such that with probability q, the coupling is suc- 
cesfull within M steps, uniformly in the initial configuration. This implies 
exponential convergence. 



4 Zhang's sandpile in infinite volume 

4.1 Definitions and main results 

In this section we work in general dimension d. We let X = [0, oo) zd de- 
note the set of infinite-volume height configurations in dimension d and 
Q = [0, l) z the set of all stable configurations. For x e Z d , the (Zhang) 
toppling operator T x is defined as 

{0 if x = y, 

Vy + hVx ti\y-x\ = l, 
rj y otherwise. 

The infinite-volume version of Zhang's sandpile model is quite different 
from its abelian sandpile counterpart. Indeed, in the infinite- volume abelian 
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sandpile model, it is shown that if a configuration can reach a stable one via 
some order of topplings, it will reach a stable one by every order of topplings 
and the final configuration as well as the topplings numbers per site are 
always the same, see [HI flOl 17]. 

In Zhang's sandpile model in infinite volume, the situation is not nearly 
as nice. Not only does the final stable realization depend on the order of 
topplings, the very stabilizability itself also does. We illustrate this with 
some examples. 

Consider the initial configuration (in d = 1) 

?? = (.. .,0,0, 1.4, 1.2,0,0,...). 

We can reach a stable configuration in any order of the topplings, but the 
final configuration as well as the number of topplings per site depend on 
which unstable site we topple first. We can choose to start toppling at the 
left or right unstable site, or to topple parallel at the same time; the different 
results are presented in Table CD 



start at site 


toppling numbers 


final configuration 


left 


(...,0,0,1,1,0,0,...) 


(...,0,0.7,0.95,0,0.95,0,...) 


right 


(...,0,1,2,3,1,0,...) 


(...,0.5,0.5,0.525,0,0.525,0.55...) 


parallel 


(...,0,0,1,1,0,0,...) 


(...,0,0.7,0.6,0.7,0.6,0,...) 



Table 1: The three possible stabilizations of (. . . , 0, 0, 1.4, 1.2, 0, 0, . . .) 



For a second example, let 

£ = (..., 0.9, 0.9, 0, 1.4, 1.2, 0, 0.9, 0.9, . . .). 

This is a configuration that evolves to a stable configuration in some orders 
of topplings, but not by others. Indeed, if we start to topple the left unstable 
site first, we obtain the stable configuration 

(. . . , 0.9, 0.9, 0.7, 0.95, 0, 0.95, 0.9, 0.9, . . .), 

but if we topple the right unstable site first, after two topplings, we reach 

£'=(..., 0.9, 0.9, 1, 0, 1, 0.6, 0.9, 0.9, . . .). 

By arguing as in the proof of the forthcoming Theorem 14.41) , one can see that 
this configuration cannot evolve to a stable configuration. 

In view of these examples, we have to be more precise about the way we 
perform topplings. In the present paper, we will use the Markov toppling 



14 



process: to each site we associate an independent rate 1 Poisson process. 
When the Poisson clock 'rings' at site x and x is unstable at that moment, 
we perform a Zhang-toppling at that site. If x is stable, we do nothing. By 
f]{t), we denote the random configuration at time t. We denote by M(x, t, rj) 
the (random) number of topplings at site x up to and including time t. 

Definition 4.1. A configuration r] G X is said to be stabilizable if for every 

x e Z d , 

lim M(x, t, rj) < oo 

a.s. In that case we denote the limit random variable by M(x, 00,77). 

Denote the collection of all stabilizable configurations by S. It is not hard 
to see that S is shift- invariant and measurable with respect to the usual Borel 
sigma field. Hence, if p is an ergodic stationary probability measure on X, 
p(S) is either or 1. 

Definition 4.2. A probability measure p on X is called stabilizable if p(S) = 
1. 

Here is our main result. 

Theorem 4.3. Let p be an ergodic translation-invariant probability measure 
on X with E M (r7o) = p < 00. Then 

1- If P > 1; then p is not stabilizable, that is, p(S) = 0. 

2. If0<p<l, then p is stabilizable, that is, p(S) = 1. 

The situation when \ < p < 1 is not nearly as elegant. Clearly, when 
we take p p to be the point mass at the configuration with constant height 
p, then p p is stabilizable for all \ < p < 1. On the other hand, we have the 
following theorem, showing that there are measures p with density close to 
I and close to 1 which are not stabilizable. 

Theorem 4.4. For all 1/2 < p < d/(2d - 1) and (2d - I)/ (2d) < p < 1, 
there exists an ergodic measure p p with density p which is not stabilizable. 

Remark There is no obvious monotonicity in the density as far as stabiliz- 
ability is concerned. Hence we cannot conclude from the previous theorem 
that for all 1/2 < p < 1 there exists an ergodic measure p p which is not 
stabilizable. 
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5 Proofs for the infinite-volume sandpile 

Starting from configuration rj = (f] x ), we define the total amount of mass 
that is lost from site x via topplings, until and including time t, by L(x, t, rj). 
We have the following identity. 

Vx{t) = Vx-L{x,t,r]) + ^ ^ L (y,t,v)- (5-1) 

\y-x\=l 

The following lemma is easy to understand, and we give it without proof. 
We write 77(00) for the (random) limiting configuration when we start with 
77. 

Lemma 5.1. For a stabilizable 77 e X , we have 

1. For every x G Z d 7 there exists a random variable L(x, 00,77) < 00 such 
that a.s. 

L(x,t,r]) t L(x, 00,7/), 

as t — > 00. 

2. The following identity holds: 

77(00) = 77(0) + AL(-, 00, 77), 
with A the matrix defined by 

{-1 ifx = y, 
h. */l2/-z| = l, 
otherwise. 

For an initial measure /x, E M and P M denote the corresponding expectation 
and probability measures in the stabilization process. We first show that no 
mass is lost in the toppling process. 

Proposition 5.2. Let /i be an ergodic shift-invariant probability measure on 
X with 

E ju(?7o) = p < 00 

which evolves according to the Markov toppling process. Then we have 

1. forO<t< oo, E M 77 (t) = p, 

2. if ji is stabilizable, then E M 7/o(oo) = p. 
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Proof. We prove 1. via the well known mass transport principle. Let the 
initial configuration be denoted by rj. Imagine that at time t = we have a 
certain amount of mass at each site, and we colour all mass white, except the 
mass at a special site x which we colour black. Whenever a site topples, we 
further imagine that the black and white mass present at that site, are both 
equally distributed among the neighbours. So, for instance, when x topples 
for the first time, all its neighbours receive a fraction l/(2d) of the original 
black mass at x, plus possibly some white mass. We denote by B(y,t) the 
total amount of black mass at site y at time t. First, we argue that at any 
finite time t, 



that is, no mass is lost at any finite time t. Indeed, had this not been the 
case, then we define t* to be the infimum over those times t for which (15.2jl 
is not true. Since mass must be subsequently passed from one site to the 
next, this implies that there is a path (x = x ,Xi, . . .) of neighbouring sites 
to infinity, starting at x such that the sites x« all topple before time £*, in 
the order given. Moreover, since t* is the infimum, the toppling times t\ of 
Xi satisfies linij^ootj = t*. Hence, for any e > 0, we can find zo so large 
that for all i > io, U > t* — e. Call a site open of its Poisson clock 'rings' 
in the time interval (t* — e,t*), and closed otherwise. This constitutes an 
independent percolation process, and if e is sufficiently small, the open sites 
do not percolate. Hence a path as above cannot exist, and we have reached 
a contradiction. It follows that no mass is lost at any finite time t, and we 
can now proceed to the routine proof of 1. via mass-transport. 

We denote by X(x, y, t, rj) the amount of mass at y at time t which started 
at x. From mass preservation, we have 




(5.2) 




(5.3) 



and 




(5.4) 



Since all terms are non-negative and by symmetry, this gives 



V?o(0) = 1^(0,3/,*, 77) 



= ^E M X(t/,0,t,r / ) = E^ (t). 



y ez d 
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To prove 2., we argue as follows. From 1. we have that for every t < oo, 
E^Q^t) = P- Using Fatou's lemma we obtain 

Mvo(oo)) = E M (lim % (t)) < liminf E„ = p, (5.5) 

and therefore it remains to show that E M (?7o(oo)) > P- We do this by letting 
a random walk run, in the same spirit as in [7]. 

Let X n , n — 0, 1, ... be a simple random walk starting at the origin. We 
write E rw , P ru , for the expectation and probability with respect to this random 
walk. Let 77 be stabilizable, and let 77 stabilize to a (random realisation) 
77(00). The numbers L(x,oo,r]) are then also fixed. For all sites x we have 
the identity 



7^(00) =7](x) -L(x, 00,77) + ^ L(y,oo,T)). 

\y-x\=l 

Plugging m x = Xk, taking expectations and observing that 



E 



I 2^ L(y,oo,r])\ =E rw (L(X k+1 ,oo,r])), 

\ \y-x k \=i j 



we obtain 

E rw {Vx k {oo) -r]x k ) = E rw (L(X k+1 ,oo,r])) - E rw (L(X k , 00, 77)). 
Now sum from k = to k = n — 1 and divide by 77 to obtain 
1 / n_1 \ 1 

-E rw y^ T] Xk (oo) -r] Xk = —E rw (L(X n , 00, 77) - L(0, 00,77)) . 

77 \ / 77 

Since L(0, 00, 77) is a fixed number, letting 77 — > 00 gives 

j n— 1 ^ n— 1 

liminf -E ri „(V'77x fe (00)) > limsup -E rw ( V ry x J (5.6) 

n— oo n f—' n ^oo 77 f—' 

K = fe=0 

Now, we choose 77 according to \i. Since /_t is ergodic, the law of 77(00) is also 
ergodic. It is well known that the scenery processes {r]x n } and {f]x n (oo)} are 
then also ergodic. Hence, the liminf and limsup are in fact limits; the right 
hand of (15.61) is equal to p, and the left hand is equal to E M (770 (00)). □ 
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Proposition 5.3. Let r](t) be obtained by the Markov toppling process start- 
ing from i] G X . Let A be a finite subset of Z rf , such that all sites in A 
toppled at least once before time t. Let (3\ be the number of internal bonds of 
A. Then 



Proof. Let (x, y) be an internal bond of A. By assumption, both x and y 
topple before time t. Suppose that x is the last to topple among x and y. As 
a result of this toppling, at least mass 1 / (2c?) is transferred from x to y and 
this mass will stay at y until time t since y does not topple anymore before 
time t. In this way, we associate with each internal bond, an amount of mass 
of at least l/(2d), which is present in A at time t. Hence the total amount of 
mass in A at time t is at least 1 / (2d) times the number of internal bonds. □ 

We can now prove Theorem 14.31 

Proof of Theorem \4-3[ We prove 1. first. Let p be any ergodic shift-invariant 
measure with E^^o) = P > 1 and suppose p is stabilizable. According to 
Proposition I5.2[ we have 



which contradicts the assumption that 77(00) is stable. 

For 2., let p be any ergodic shift-invariant probability measure on X with 
E M (r/o) = P < 21 an d suppose that p is not stabilizable. We will now show 
that this leads to a contradiction. 

Define C n (t) to be the event that before time t, every site in the box 
[— n,n] d topples at least once. Since p is not stabilizable we have that 
P M (C n (£)) —>■ 1 as t —>■ 00. Indeed, if a configuration is not stabilizable, 
all sites will topple infinitely many times as can be easily seen. 

Choose e > such that 1 — e > 2p. For this e, there exists a non-random 
time T e > such that for all t > T e , 



From Proposition 15.31 we have that at time t > T € > 0, with probability at 
least 1 — e, the following inequality holds: 




(5.7) 



E m Moo)) = E m (t7 ) = p>1, 



(5.8) 



P M (C n (t)) >l-e. 



(5.9) 



E x e[-n,n] d7 l*( t ) . 1 (2n) d 



(5.10) 



(2n + l) d ~ 2(2n + l) d ' 
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Therefore, we also have 



^ ' Es e[ -n,n]^(*A > 1 (1 (2n) 



(2n + l) d / ~ 2 V ; (2n + 
Since 2p < 1, we can choose n so large that 

v (2n) d 

Using the shift-invariance of p and the toppling process, for t > T e , we find 

„ m „ /E l6[ - n ,f1x(i)\ I (2n) d 

E u 7?n(t) = E u 7 i — -^-i — > -(1 - e)-^ — —-: > p. (5.11) 

M/uw M (2n + l) d y ~ 2 V ; (2n + l) d F y ' 

However, from Proposition 15.21 we have for any finite t that E At r/o(t) = 
E^ (0) =p. □ 

Proof of Theorem \4.4\ We start with p > (2d — l)/(2d). To understand the 



idea of the argument, it is useful to first assume that we have an unstable 
configuration rj on a bounded domain A (with periodic boundary conditions) 
with the property that rj x > 1 — l/(2rf), for all igA. On such a bounded 
domain, we can order the topplings according to the time at which they 
occur. Hence we can find a sequence of sites Xi,x 2 , ■ ■ ■ (not necessarily all 
distinct) and a sequence of times t± < t 2 < ■ ■ ■ such that the i-th toppling 
takes place at site ij 6 A at time t$. At time ti, x\ topples, so all neighbours 
of x\ receive at least l/(2d) from x\. This means that all neighbours of 
x\ become unstable at time t±, and therefore they will all topple at some 
moment in the future. As a result, x\ itself will also again be unstable after 
all its neighbours have toppled, and hence x\ will topple again in the future. 

In an inductive fashion, assume that after the k-th toppling (at site Xk 
at time tk), we have that it is certain that all sites that have toppled so far, 
will topple again in the future, that is, after time tk- Now consider the next 
toppling, at site x k+ i at time t k+ i. If none of the neighbours of x k +i have 
toppled before, then a similar argument as for X\ tells us that x k +i will topple 
again in the future. If one or more neighbours of x k+ i have toppled before, 
then the inductive hypothesis implies that they will topple again after time 
tfc+i- Hence, we conclude that all neighbours of x k +i will topple again which 
implies, just as before, that x k +i itself will topple again. We conclude that 
each sites which topples, will topple again in the future, and therefore the 
configuration can not be stabilized. 

This argument used the fact that we work on a bounded domain, since 
only then is there a well-defined sequence of consecutive topplings. But with 
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some extra work, we can make a similar argument work for the infinite- 
volume model as well, as follows. 

Let so > be so small that the probability that the Poisson clock at 
the origin 'rings' before time sq is smaller than the critical probability for 
independent site percolation on the <i-dimensional integer lattice. Call a 
site open if its Poisson clock rings before time s . By the choice of s , all 
components of connected open sites are finite. For each such component of 
open sites, we now order the topplings that took place between time and 
time So- For each of these components, we can argue as in the first paragraphs 
of this proof, and we conclude that all sites that toppled before time sq, will 
topple again at some time larger than s . We then repeat this procedure 
for the time interval [s ,2so]; [2so,3s ], • • •> an d conclude that at any time, 
a site that topples, will topple again in the future. This means that the 
configuration is not stabilizable. Hence, if we take a measure p p such that 
with /ip-probability 1, all configurations have the properties we started out 
with, then p p is not stabilizable. 

Next, we consider the case where 1/2 < p < d/(2d — 1). Consider a 
measure p p whose realizations are a.s. 'checkerboard' patterns in the following 
way: any realization is a translation of the configuration in which all sites 
whose sum of the coordinates is even obtain mass 2p, and all sites whose sum 
of coordinates is odd obtain zero mass. Consider a site x with zero mass. 
Since all neighbours of x are unstable, these neighbours will all topple at some 
point. By our choice of p, x will become unstable precisely at the moment 
that the last neighbour topples - this follows from a simple computation. 
By an argument pretty much the same as in the first case, we now see that 
all sites that originally obtained mass 2p, have the property that after they 
have toppled, all their neighbours will topple again in the future, making 
the site unstable again. This will go on forever, and we conclude that the 
configuration is not stabilizable. Hence, p p is not stabilizable. □ 

Remark The argument in case of parallel topplings is simpler and works for 
all p > 1/2; in that case the checkerboard pattern is preserved at all times, 
preventing stabilization. 
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